Muon pair creation by two x-ray laser photons in the field of an atomic nucleus 
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The generation of muon-antimuon pairs is calculated in the collision of an ultrarelativistic bare 
ion with an intense x-ray laser beam. The reaction proceeds nonlinear ly via absorption of two laser 
photons. By systematic study throughout the nuclear chart we show that the interplay between the 
nuclear charge and size, along with the possibility of nuclear excitation leads to saturation of the 
total production rates for high-Z ions, in contrast to the usual Z 2 scaling for pointlike projectiles. 
The process is experimentally accessible by combining present-day ion accelerators with near-future 
laser sources and in principle allows for the measurement of nuclear form factors. 

PACS numbers: 12.20.Ds, 21.10.-k, 32.80.Wr, 34.90.+q 
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In recent years the interest in electron-positron pair cre- 
ation in combined laser and nuclear Coulomb fields has 
been stimulated by the large advances in high-power laser 
devices and ion accelerators (see [3, 0, H| and references 
therein). The process is interesting since it proceeds by 
the absorption of real photons from the laser field and 
a virtual photon from the Coulomb field. Most stud- 
ies have considered the low-frequency tunneling regime, 
where laser field strengths close to the Schwinger value 
E^ = m 2 c 3 /he = 1.3 x 10 16 V/cm are required. At this 
field intensity, the electric work performed on an elec- 
tron of charge e along the Compton wavelength = 
h/m e c equals its rest energy m e c 2 . In the opposite high- 
frequency multiphoton regime, pairs are created via the 
simultaneous absorption of several high-energy photons 
with huj ~ m e c 2 from a rather weak laser field. In both 
cases, the presence of the nucleus is required to guarantee 
energy-momentum conservation. 

The field strengths and frequencies required for e + e~ 
pair creation in the respective regimes are four orders 
of magnitude larger than provided by the most-advanced 
present laser systems. This gap can be bridged in laser- 
ion collisions, where the ions are counterpropagating the 
laser beam at relativistic speed [H, 0, [1] . Then the laser 
field strength and frequency as seen by the nucleus in its 
rest frame are enhanced by the relativistic Doppler shift. 
The only experimental observation of e + e~ pair produc- 
tion in laser and Coulomb fields has been achieved at 
SLAC (Stanford, California) in collisions of ultrarelativis- 
tic electrons with an intense optical laser pulse The 
pairs were generated indirectly via Compton backscat- 
tering and a subsequent collision of real photons. 
In view of the ongoing technical progress the question 
arises, whether not only e + e~ but also fj, + p~ pairs can 
be produced with the emerging near-future laser sources. 
Muon production in the tunneling regime appears rather 
hopeless, though, since the required field needs to be 
close to E ( / } = p 2 E ( g ) = 5.6 x 10 20 V/cm, with the 
muon-to-electron mass ratio p = m^/m e f=a 207. Even by 
exploiting the 7-factors which will soon be provided by 
the Large-Hadron Collider at CERN (Geneva, Switzer- 



land) this value seems out of reach (7 = 7000 for protons 
[Bj]). As an alternative, it has recently been proposed 
that p~ pairs could be produced in laser and Coulomb 
fields within a two-step process Q: First, an e + e~ pair 
is created via tunneling in a laser-ion collision, and after- 
wards driven by the laser field into an annihilating rec- 
ollision e + e~ — > The latter process has also been 

studied in 0] , with the electron and positron originating 
from a positronium atom or an e + e~ plasma. 
In this Letter, we calculate direct p + \i~ production in 
highly energetic laser-ion collisions. Rather than the tun- 
neling regime, we consider multiphoton muon production 
via absorption of few high-frequency laser photons. The 
latter are assumed to stem from an x-ray free electron 
laser (XFEL) beam with fkj — 12keV [7|, which col- 
lides head-on with an ultrarelativistic nucleus moving at 
7 = 7000. In the nuclear rest frame the photon energy 
amounts to tkv « 2jhu>o = 168 MeV, so that the energy 
threshold Ae = 2m At c 2 for production can be over- 

come by two-photon absorption from the XFEL field 0] . 
Corresponding large-scale XFEL facilities are presently 
being developed at SLAC and DESY (Hamburg, Ger- 
many), where peak intensities close to 10 20 W/cm 2 are 
envisaged 0- The projectile nucleus is modelled by an 
extended charge distribution, whose shape is shown to 
have significant impact on the muon creation rates. 
At first sight, e + e~ and p + pT pair production in com- 
bined laser and Coulomb fields seem to be very simi- 
lar processes since the electron and muon only differ by 
their mass (and lifetime). In this picture, the corre- 
sponding production probabilities would coincide when 
the laser field strength and frequency are scaled in ac- 
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cordance with the mass ratio, i.e. W. 
W e+e -{E^ e \Lu^) for EM = p 2 E^ and uW = pu 
This simple scaling argument does not apply, however, 
as the large muon mass is connected with a correspond- 
ingly small Compton wavelength "Xq = %q / p ~ 1.86 fm, 
which is smaller than the radius of most nuclei. As a re- 
sult, the nucleus does not look pointlike to the muon and 
its finite extension must be taken into account. Pro- 
nounced nuclear size effects have also been found for 
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production by single 7-photon impact on nuclei 
and in relativistic heavy- ion collisions [ToL 11 1 . 
Following the usual theoretical approach to lepton pair 
creation in combined laser and Coulomb fields (see, e.g., 
0), we write the process amplitude in the nuclear rest 
frame as 
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The muons are created with free momenta p+ and de- 
scribed by relativistic Volkov states \& p± [12|] which in- 
clude their interaction with the laser field to all orders. 
The nuclear field V(r) is taken into account within the 
first-order of perturbation theory. We assume the nucleus 
to be spherically symmetric with a Gaussian charge dis- 
tribution 
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where Z is the atomic number and the parameter a is re- 
lated to the nuclear rms charge radius by r rms = y/S]2a. 
The charge density ^ will serve us for a systematic study 
of the nuclear size effect on the muon creation process. 
We note that light nuclei with mass number A ~ 10 are 
Gaussian-shaped to a good approximation, whereas the 
charge density of heavy isotopes (A ~ 100) is usually 
parametrized by a Fermi distribution [l3j ]. The electro- 
static potential generated by the charge density ((2|) is 
V(r) — (Ze/r)$(r/a), with the error function The 
amplitude in Eq. (TT]) can be evaluated analytically by ex- 
panding its periodic part into a Fourier series. This way, 
one arrives at the integral 



V(q) = J d 3 r V(r) e^ r = J e ~ii-/^) 2 



(3) 



where q = p + + p — nk is the momentum transfer to 
the nucleus. Within the multiphoton regime, the laser- 
dressed momenta p± introduced here practically coincide 
with the free momenta p± [l2j. Moreover, n = 2 laser 
photons of momentum k are absorbed from the field. 
The first factor on the right-hand side of Eq. ([3]) is the 
Fourier transform of the Coulomb potential of a point- 
like nucleus. The second factor is the elastic nuclear form 
factor F(q 2 ), which describes the correction due to the fi- 
nite nuclear size and leads to substantial reduction of the 
process probability when q ~ m^c approaches or exceeds 
the inverse nuclear radius. The fully differential rate for 
two-photon muon pair creation is obtained by summing 
the square of the amplitude ([1]) over the lepton spins: 



adopts the structure dR c \ = dR Z 2 F 2 (q 2 ), with dRo 
denoting the differential production rate for a pointlike 
proton. This structure is analogous, e.g., to the form 
factor-corrected Rutherford cross section or the Rosen- 
bluth formula for electron-nucleon scattering, which also 
factorize into a product of the result for a point scat- 
terer and the form factor squared [lj|. An additional 
contribution to creation comes from the inelastic 

channel, where the nucleus is excited due to the recoil 
imparted on it [U [3] • The Z protons inside the nucleus 
act incoherently here, and one has to a good approxima- 
tion dRind w dR Z(l — F 2 (q 2 )). The inelastic process 
becomes important when f 2 C 1. In principle, both pro- 
duction channels are separately accessible in experiment 
by observing the final nuclear state in coincidence. 
Figure 1 demonstrates the nuclear size effect by showing 
angle-differential rates for elastic muon pair creation in 
ultrarelativistic collisions of various ions with an intense 
XFEL beam of linear polarization. The emission angle 
refers to the laser wave vector and the rates are scaled 
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by 1/Z 2 . The laser intensity parameter £ M = n 
lates to the muon mass; the value chosen corresponds to 
2.54 x 10 22 W/cm 2 . We point out that this intensity ex- 
ceeds the original design values at SLAC and DESY by 
2-3 orders of magnitude, but can be attained by antic- 
ipated facility extensions and improved x-ray focussing 
techniques [15j . The respective rms charge radii are 
0.875 fm (proton), 2.470 fm ( 12 C), 3.929fm ( 64 Zn), and 
5.851 fm ( 238 U) [16[. For comparison we also show the 
result for a pointlike nucleus. With increasing nuclear 
size, the production rates are strongly reduced, and the 
maximum of the angular distributions is shifted towards 
smaller angles. The latter is because the nuclear form fac- 
tor cuts the contributions from large momentum trans- 
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Equation accounts for the elastic channel of the pro- 
cess, where the nucleus remains in its ground state. It 



FIG. 1: Angular spectra of one of the muons produced by 
two-photon absorption from an intense XFEL beam (huo = 
12keV, ^ = 6.8 x 1(T 5 ) colliding with various ultrarelativistic 
nuclei (7 = 7000), as indicated. The rates are taken in the 
nuclear rest frame and scaled by the nuclear charge. 
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fers, which otherwise give rise to large emission angles. 
The reduction of the total rates with respect to the result 
for a point charge arises from the fact that the particles 



are created at typical distances 



When this dis- 



tance is smaller than the nuclear radius, only a fraction 
eZ e g = 4-7T J r 2 g(r)dr of the total nuclear charge con- 
tributes effectively to the process. 

Figure 2 shows integrated rates of production for 

several projectiles. The nuclear isotope of largest abun- 
dance was chosen always. The values again refer to the 
nuclear rest frame. For a point proton, the total rate 
would amount to i?o = [oi 2 /IQ^^loT^) « 0.58s _1 
(cp. Eq. (26) in 0), displaying the typical dependence 
of a two-photon process. Here, a denotes the finestruc- 
ture constant and the dimensionless function is of 

the order of 0.1 in the relevant frequency range (i.e., 
1 < frui/m^c 2 < 2). As compared to the point-charge 
result, the elastic production rates are reduced by 23% 
for proton impact and by 77% for 12 C impact. For the 
heavier projectiles, the reduction factors amount to 0.13 
( 27 A1), 0.087 ( 40 Ca), 0.05 5 ( 64 Zn), 0.042 ( 84 Kr), 0.026 
( 120 Sn), 0.020 ( 142 Nd), 0.010 ( 208 Pb), and 7.1 x 10~ 3 
( 238 U), respectively. Since the elastic process probability 
increases with the projectile charge as Z 2 and decreases 
with its size as ~ exp[— ^(r lms /%^) 2 ], a maximum elas- 
tic rate arises for atomic numbers around Z as 60. This 
is in contrast to e + e~ pair creation which increases as 
Z 2 throughout [H, 0, Q , apart from higher-order correc- 
tions in aZ which slightly modify this behaviour at high 
Z. Note that the latter corrections are of minor impor- 
tance in the present situation, as the muons are produced 
well above the energetic threshold. The emergence of the 
maximum in Fig. 2 can be understood by consideration of 
the effective nuclear charge Z r contained in a sphere of 
radius . The corresponding values 3.66 ( 84 Kr), 3.83 
( 120 Sn), 3.96 ( 142 Nd), 3.96 ( 208 Pb), and 3.73 ( 238 U) agree 
with the observed location of the maximum. This simple 
picture also explains a local minimum for atomic numbers 
around Z « 70 (not shown), since Z c g « 3.72 for 174 Yb. 
We stress that the maximum remains at the same posi- 
tion, when the nuclei are modelled by uniformly charged 
spheres, which represents a better approximation than 
Eq. ^ to the Fermi charge distribution of the relevant 
heavy isotopes where the maximum occurs 13[. A max- 
imum also arises for elastic muon production by a single 
photon of twice the energy. Figure 2 moreover displays 
the total rate i?tot = i? c i + -Rmci, which does not exhibit a 
maximum but saturates at high Z values since the contri- 
bution from the inelastic channel increases with nuclear 
charge. The total rates are still considerably smaller than 
the point-proton results: The reduction factors amount 
to, e.g., 0.35 ( 12 C), 0.13 ( 40 Ca), 0.069 ( 84 Kr), 0.036 
( 142 Nd), and 0.018 ( 238 U). For e+e~ pair creation, the 
inelastic channel is negligible since q ~ m e c <C h/r rms , so 
that F 2 w 1 and i? to t « RqZ 2 . 
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FIG. 2: Nuclear- frame rates for muon creation in the XFEL- 
nucleus collisions as in Fig. 1. The squares and triangles show 
the total and elastic rates for various nuclei, respectively, con- 
nected by fit curves. The dotted line oc Z 2 holds for a point 
nucleus. 
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FIG. 3: Angular spectra in the laboratory frame for two- 
photon muon pair production by 12 C nuclei of different shapes 
and an XFEL beam. The upper (lower) two curves show total 
(elastic) rates. The collision parameters are as in Fig. 1. 



Figure 3 shows angular spectra in the laboratory frame 
for carbon impact. Apart from the Gaussian charge dis- 
tribution of Eq. ((2|), a uniform density within a sphere 
of radius ^/5/3 r rms was assumed, which has equal rms 
radius. As compared to the result for a Gaussian-shaped 
12 C nucleus, the elastic (total) production rate is smaller 
by 13% (7%) for the hard sphere. The relative difference 
grows for heavier nuclei, reaching about 40% (20%) for 
238 U. This demonstrates that the muon production pro- 
cess is sensitive not only to the nuclear radius, but also to 
the nuclear shape. In principle, the process might even be 
used for determination of the nuclear form factor, which 
is traditionally accomplished by electron scattering. 
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According to our results, muon pair creation in XFEL- 
nucleus collisions is experimentally feasible by employing 
near-future technology. In the laboratory frame, the pro- 
duction rates of Fig. 2 are reduced by a factor 1/7 due to 
relativistic time dilation and, accordingly, reach values of 
about i?iab ral0 _2 s _1 . In the collision of an ion beam 
containing N — 10 11 particles Q with an XFEL pulse 
of t = 100 fs duration [7j, the probability for production 
of one muon pair is W„+„- ~ R\ a bTN/2 10~ 4 . At 
the envisaged average XFEL repetition rate of 40 kHz, 
one muon production event per second is obtained. The 
muon yield scales with the laser intensity squared. The 
typical lab-frame energy of the muons ~ 500 GeV is 
highly relativistic; the muonic lifetime is accordingly in- 
creased to ~ 10 ms. e + e~ pairs are also produced in the 
collision by single-photon absorption in the nuclear field 
(Bethe-Heitler process [12|). We stress, however, that the 



subsequent reaction e + e 



MM 



via laser-driven recol- 



lisions [3J does not occur in the multiphoton regime con- 
sidered here, because the relative electron-positron mo- 



mentum Ap 



satisfies -C Ap <C 2to m c. Detected 



muons have thus been produced via the direct process 
A + riLj — > A + fi + pT , with n = 2. Higher photon orders 
(n > 3) are suppressed by an additional factor ^ < 1. 
In XFEL-nucleus collisions, also bound-free p + M~ P arr 
production can take place, with the fi~~ being created in 
a bound state of the atomic nucleus. In order to esti- 
mate the corresponding rate, we apply an appropriatly 
modified version of the theory developed in |l7| for e + e~ 
production with K-shell capture. Assuming the same col- 
lision parameters as in Fig. 2 and a pointlike proton, we 
obtain a rest-frame rate of 1.94 x 10~ 6 s^ 1 . This value is 
by five orders of magnitude smaller than the correspond- 
ing rate for free muon pair creation in Fig. 2. For point 
projectiles, the bound-free production channel raises like 
Z 5 with atomic number [17|. This scaling is also consid- 
erably damped by the finite nuclear size. 
Muon pair creation has been studied before in relativis- 
tic heavy-ion collisions [10]. The process is formally re- 
lated to muon production by a single photon in the field 
of an atomic nucleus 0, [l8| by the Weizsacker- Williams 
method, which describes the transverse electromagnetic 
field of an ultrarelativistic nucleus by a spectrum of vir- 
tual photons [HI]. In accordance with the present re- 
sults, the finite nuclear extent was found to reduce the 
production probability substantially. In the high-energy 
limit, the argument of the famous Bethe-Heitler loga- 
rithm log(2fi<x> /m M c 2 ) is modified by the factor Xq' V r rms 



loL This factor is also responsible for the rate sup- 
pression of two-photon muon pair creation here. Numer- 
ical calculations of bound-free muon pair production in 
238 U- 238 U collisions have found a cross section reduction 
by five orders of magnitude due to the nuclear extent [l|| . 
It is interesting to note that exclusive dilepton produc- 
tion in hadron collisions also represents a search tool for 
physics beyond the standard model [20| . 



In conclusion, direct production of pairs by two- 

photon absorption from a high-frequency laser wave col- 
liding with an atomic nucleus was calculated. The pro- 
cess is sensitive to the nuclear form factor. It could 
be realized experimentally by combining the radiation 
from upcoming XFEL sources with an ultrarelativistic 
ion beam from the present generation of heavy-ion accel- 
erators. The highest elastic production rates are reached 
by projectiles in the lanthanoid region (Z rs 60), whereas 
the total production rates saturate at high Z values. 
Useful input by T. J. Burvenich, K. Z. Hatsagortsyan, 
and A. B. Voitkiv is gratefully acknowledged. 
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